Auxin regulates source-sink carbohydrate partitioning and reproductive organ development in rice

Significance Effective communication between source organs and sink organs is pivotal in carbohydrate assimilation and partitioning during plant growth and development. Auxin is required for many aspects of plant growth and development. However, very little is known about how these two important classes of molecules coordinate and co-regulate plant developmental processes. In this study, we elucidate an OsARF18-OsARF2-OsSUT1–mediated auxin signaling cascade regulating carbohydrate partitioning between the source and sink tissues in rice, which is essential for proper development of rice reproductive organs. Our findings represent a major step forward in increasing our knowledge of sucrose transport regulation in plants and have important implications in improving crop yield through better coordination of source and sink activities.

In higher plants, the sink organs, such as flowers, fruits, and seeds, are heterotrophic in nature and rely on nutrients supplied from the photosynthetically active organs (e.g., leaves, termed source organs) for their growth and development (1)(2)(3)(4). Higher plants utilize the phloem sieve elements for long-distance transport of nutrients (mainly sucrose) from the source to the sink organs. The turgor pressure generated through the osmotic effect of sucrose loading into the phloem at the source and the unloading at the sink creates the driving force for long-distance translocation of all other compounds, including nutrients, water, and signaling molecules in the phloem (3,5).
Sucrose translocation is mainly mediated by two sucrose transporter families: Sugars Will Eventually be Exported Transporters (SWEETs) and sucrose transporters or sucrose carriers (SUTs/SUCs). SWEET proteins act to export sucrose and hexose to the phloem apoplasm (6,7), while SUT proteins mainly import sucrose from the apoplasm into the sieve tube symplasm (8)(9)(10). SUT proteins have 12 transmembrane domains that form a pore in the membrane to allow the passage of sucrose, and they utilize the energy stored in the proton gradient across the membrane generated by H + -adenosinetriphosphatases to drive sucrose transport (3,11). In addition, cell wall invertases hydrolyze the unloaded sucrose from the phloem in the sink tissues into glucose and fructose, which are subsequently taken up by membrane-bound monosaccharide transporters of the recipient cells in developing tissues (5,12). Disruption of the various sugar transporters often results in elevated sugar content in the source leaves and abnormal reproductive development (including abnormal anther development, impaired male fertility, and seed filling) because of carbohydrate deprivation (10,(13)(14)(15)(16)(17).
Given the important roles of sugar transport in plant growth and development, it is not surprising that the activities of various sugar transporters must be tightly regulated. Emerging evidence suggests that sugar transporters can be regulated at both the transcriptional and posttranscriptional levels (phosphorylation and dephosphorylation) (18). For example, the Carbon Starved Anther (CSA) gene encodes a putative R2R3 MYB-type transcription factor regulating the expression of OsMST8 during male reproductive development (17), and the rate of carbon export from source leaves is controlled by ubiquitination and phosphorylation of SUC2 (19). There is also evidence indicating that biotic and abiotic factors (including light, water and salt stress, and pathogen attack) and phytohormones have an impact on the expression of sugar Significance Effective communication between source organs and sink organs is pivotal in carbohydrate assimilation and partitioning during plant growth and development. Auxin is required for many aspects of plant growth and development. However, very little is known about how these two important classes of molecules coordinate and co-regulate plant developmental processes. In this study, we elucidate an OsARF18-OsARF2-OsSUT1-mediated auxin signaling cascade regulating carbohydrate partitioning between the source and sink tissues in rice, which is essential for proper development of rice reproductive organs. Our findings represent a major step forward in increasing our knowledge of sucrose transport regulation in plants and have important implications in improving crop yield through better coordination of source and sink activities.
transporters (18,20). For example, several studies have suggested that auxin signaling might be linked to sugar metabolism (4,21). Down-regulation of the tomato IAA9 gene and auxin response repressor ARF4 leads to increased input of sugar to the fruit (22,23). Despite the progress made in this area, however, the molecular mechanisms regulating carbohydrate partitioning between source and sink tissues remain elusive.
In a previous study, we reported that the rice dao mutant, which is defective in a gene encoding a 2-oxoglutarate-depenedent-Fe (II) dioxygenase, failed to convert indole-3-acetic acid (IAA) into 2-oxoindole-3-acetic acid (OxIAA), the important step toward the degradation of IAA (24,25). The dao mutant displayed a pleiotropic reproductive phenotype, including closed spikelets, indehiscent anthers, and development of parthenocarpic seeds filled with a sucrose-rich liquid. In this study, we show that DAO-mediated auxin homeostasis plays an important role in regulating carbohydrate partitioning between the source and sink tissues to regulate reproductive organ development in rice.

Elevated Auxin Levels Cause Closed Spikelets, Indehiscent
Anthers, and Formation of Parthenocarpic Seeds. We previously showed that the rice dao mutant displays parthenocarpic seeds due to elevated auxin levels (24). To test whether other reproductive defects (closed spikelets and indehiscent anthers) of the rice dao mutant are also caused by elevated auxin levels in these organs, we generated transgenic rice plants overexpressing OsYUC4 (SI Appendix, Fig. S1 A and B), a member of the OsYUC (YUCCA) gene family that is highly expressed in rice anthers (26). The YUC family of flavin monooxygenases catalyzes the rate-limiting step in IAA biosynthesis (27). As expected, the OsYUC4-overexpressing plants exhibited a daolike phenotype, including closed spikelets, indehiscent anthers, and parthenocarpic seeds. In addition, treatment of wild-type (WT) panicles with exogenous auxin (2,4-D) also led to similar phenotypes ( Fig. 1 A-D and SI Appendix, Fig. S1 C-J). These results indicate that elevated auxin levels cause abnormal development of reproductive organs in rice.
The Lodicules Fail to Enlarge at Anthesis in the dao Mutant.
The closed spikelets phenotype of the dao mutant, OsYUC4 overexpressors, or WT plants treated with exogenous auxin led us to analyze lodicule development, which is a grass-specific floral organ with scale-like shapes that swell rapidly immediately before anthesis to promote spikelet opening, allowing subsequent pollination/fertilization (28)(29)(30). No significant difference was observed in the sizes of WT and dao mutant lodicules before flowering (Fig. 1 E and F). At anthesis, the size of lodicules in the WT was obviously larger than that in the dao mutant (Fig. 1G). After pollination, the lodicules in the WT became shrunken, leading to closed glumes. In contrast, the lodicules failed to expand during anthesis, and the glumes remained closed in the dao mutant (Fig. 1G).
Abnormal Sugar Partition in the dao Mutant. The failure of lodicules to expand at anthesis and the closed glume phenotype of the dao mutant prompted us to speculate that this phenotype might be caused by abnormal cell wall expansion of the lodicules, defects in water flux, and/or carbohydrate partitioning. Thus, we first analyzed the lodicule morphology of WT and dao mutants by scanning electron microscopy. No apparent differences were found in the cell wall of lodicules in WT and the dao mutant (SI Appendix, Fig. S2). As concentration of sucrose in the phloem is thought to be the dominant osmoticum that drives translocation of all other solutes, including water flux (3, 7), we measured the sugar content in the lodicules of WT and dao mutants. As expected, sucrose levels were significantly lower in the lodicules of dao compared to those in WT (Fig. 1H). In addition, we also measured the sugar contents in the sink tissues (anthers and ovaries), source tissues (flag leaves), and surrounding tissues (stems) of both WT and dao mutants. Strikingly, sucrose levels were significantly lower in the sink tissues of dao compared to those in WT. In contrast, higher levels of sucrose were found in the source tissues of the dao mutants compared to those in WT ( Fig. 1 I-L). These observations indicate a potential defect in carbohydrate partitioning between the source and sink tissues in the dao mutant.
To directly test this possibility, we performed isotopelabeling experiments by feeding the flag leaves of WT and dao mutants with [ 14 C] sucrose (Fig. 1M). Radioactivity testing showed that 24 h after treatment, the dao mutants had more labeled sucrose in the flag leaves and stems than the WT. In contrast, the WT had more [ 14 C] label in the anthers and ovaries than the dao mutant plants ( Fig. 1 N and O). Additionally, we fed the flag leaves of WT and dao mutants with the same concentration of 5, 6-carboxyfluorescein diacetate (5, 6-CFDA), a membrane-permeable, nonfluorescent dye (SI Appendix, Fig.  S3A). Upon entering the cell, 5, 6-CFDA is converted into a fluorescent, membrane-impermeable fluorescent form (CF) tracer. In the symplasm, it can be translocated through the phloem (31). Half an hour after the treatment, strong fluorescence was detected in the sink tissues (panicles, spikelets, and lodicules) of WT but not of the dao mutant plants ( Fig. 1 P-S and SI Appendix, Fig. S3 B-D). In addition, we submerged the excised stem of WT and dao mutant in water containing 5, 6-CFDA. After 30 min of treatment, the CF tracer was detected in the stem of both the WT and dao mutant (SI Appendix, Fig. S4 A-F). Moreover, we found that 5, 6-CFDA movement in the panicles and lodicules of WT plants treated with 2,4-D (auxin applied to the flag leaf) was inhibited by auxin treatment (SI Appendix, Fig. S4 G-L). Together, these results suggest that defects in sugar loading in the dao mutant are mainly caused by elevated auxin levels rather than defective phloem elements.

Mutations in Ossut1
Cause Defective Development of the Lodicule, Anther, and Ovary. To investigate the molecular mechanisms of auxin-regulating source-sink carbohydrate partitioning in rice, we conducted RNA sequencing (RNA-seq) analysis of the lodicules in the WT and dao mutant (lodicules collected at 9:00 to 10:00 AM before flowering). A total of 1,144 differentially expressed genes (DEGs) were detected, among which 174 genes were up-regulated and 970 were down-regulated ( Fig. 2A). Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway ( Fig. 2B and Dataset S1) and gene ontology analysis indicated enrichment of plant hormone signal transduction, carbon metabolism, starch, and sucrose metabolism in the DEGs (SI Appendix, Fig. S5 and Dataset S2). Notably, Os03g0170900 (known as OsSUT1) was among the 10 most significantly down-regulated DEGs in the dao mutant ( To test whether the altered expression of OsSUT1 might contribute to the observed phenotype of the dao mutant, we generated OsSUT1 RNA interference knockdown transgenic plants (SI Appendix, Fig. S6D). These transgenic plants showed  the spikelet closing and anther indehiscence phenotype (SI Appendix, Fig. S6 E and F). In addition, we obtained Ossut1 knockout mutants using the CRISPR-Cas9 technique (Fig. 2G). As expected, the Ossut1 knockout mutants exhibited a dao-like phenotype including closed spikelets, defective pollen germination, and formation of parthenocarpic seeds ( Fig. 2H and SI Appendix, Fig. S6 G-N). Moreover, measurements of sugar content using gas chromatography-mass spectrometry (GC-MS) showed that the Ossut1 mutant had significantly lower levels of sucrose in the sink tissues but significantly higher levels of sucrose in the source and surrounding tissues compared with the WT (Fig. 2 I-L). Feeding experiments with [ 14 C] sucrose revealed results similar to those of GC-MS experiments ( Fig. 2 M and N). These results demonstrate that OsSUT1 plays a critical role in carbohydrate partitioning during rice reproductive development.
Knocking Out OsARF2 Phenocopies the dao Mutant. Auxin is a major plant hormone that regulates many aspects of plant growth and development, including embryogenesis, seed development, root development, seedling growth, vascular patterning, and reproductive development (26,(32)(33)(34). The canonical auxin-signaling pathway is centered on activation or repression of gene expression by a family of auxin response factors (ARFs) (35,36). The rice genome encodes 25 OsARF proteins (37).
To identify the OsARFs that may contribute to the dao phenotypes, we analyzed the RNA-seq data of the 25 OsARFs in the lodicules of WT and the dao mutant (SI Appendix, Fig. S7A). We found that expression of OsARF2 was significantly lower in the dao mutants ( Fig. 3A and SI Appendix, Fig. S7B), suggesting that expression of OsARF2 might be negatively regulated by the elevated level of auxin in the dao mutant. Consistent with this notion, OsARF2 expression was significantly decreased in the OsYUC4 overexpression lines and in WT panicles treated with exogenous auxin (Fig. 3 B and C, and SI Appendix, Fig.  S7 C and D). qRT-PCR experiments demonstrated that OsARF2 was mainly expressed in stems, young anthers, and developing ovaries (SI Appendix, Fig. S7E).
To test whether OsARF2 plays a role in regulating reproductive development, we performed targeted mutagenesis of OsARF2 in the ZH11 cultivar background using CRISPR-Cas9 technology. We obtained two independent, Cas9-free, homozygous Osarf2 mutants, Osarf2-1 and Osarf2-2 (Fig. 3D). No significant differences were observed between the Osarf2-1 mutant and WT plants during vegetative growth. As in the dao mutant, some glumes of the Osarf2-1 mutants failed to open normally at the anthesis stage (Fig. 3 E-G). In addition, some anthers remained completely indehiscent and did not release mature pollen grains, and some seeds were parthenocarpic (SI Appendix, Fig. S7 F-H).
To further confirm that these phenotypes were caused by mutations in the OsARF2 gene, we performed a functional complementation experiment. We introduced a 6.8-kb WT DNA fragment containing the OsARF2 promoter and its coding region into the Osarf2 mutant. The transgene fully rescued the Osarf2 mutant phenotypes, including spikelet opening, anther dehiscence, and seed development ( Fig. 3H and SI Appendix, Fig. S7I). Collectively, these results indicate that OsARF2 plays an essential role in regulating spikelet opening, anther dehiscence, and seed development.

OsARF2 Directly Regulates the Expression of OsSUT1 via
Binding to the Sugar-Responsive Elements (SuREs). ARF proteins are known to regulate downstream gene expression by binding to the cis-element motifs termed auxin-responsive elements (AuxREs: TGTCNN) in the target gene promoters using a conserved DNA-binding domain at the N terminus (38)(39)(40)(41). We identified eight AuxREs and four sugar-responsive elements (SuRE1-SuRE4: GTCTC) (42) in the OsSUT1 promoter region (Fig. 3 I and J). Interestingly, the sequences of AuxREs and SuREs are highly similar. To test whether OsARF2 directly regulates OsSUT1 expression, we prepared various deletion fragments of the OsSUT1 promoter (from À2,700 bp to À84 bp), fused them to the luciferase reporter gene coding sequence, and performed transient expression experiments in rice protoplasts. LUC activity was significantly up-regulated to a comparable level when OsARF2 was coexpressed with the P-S1::LUC to P-S7::LUC reporters, but not with the P-S8::LUC reporter (Fig. 3I), suggesting that OsARF2 may regulate OsSUT1 expression by direct binding to the SuRE1 or SuRE2 element, both of which are present in the S1-S7 fragments but missing in the S8 fragment. Chromatin immunoprecipitation (ChIP)-qPCR assays showed that, indeed, OsARF2 could directly bind to the CH1 and CH2 fragments that contain the SuRE1 or SuRE2 element, respectively, but not to other fragments (Fig. 3J). Moreover, electrophoretic mobility shift assay (EMSA) also showed that OsARF2 could directly bind to the SuRE1 and SuRE2 fragments (Fig. 3K), and the binding activity was gradually reduced by increased concentrations of unlabeled probes (Fig. 3 L and M). Further, mutation or deletion of the SuRE1 and SuRE2 elements reduced or abolished the binding activity of OsARF2 (SI Appendix, Fig. S8 A-C). These results verify that OsARF2 directly regulates OsSUT1 expression through direct binding to the SuRE1 and SuRE2 elements in the OsSUT1 promoter.
The observation that OsARF2 directly binds to the SuRE1 and SuRE2 elements but not the AuxRE elements in the OsSUT1 promoter (Fig. 3K) suggested to us that SuREs might represent a novel type of ARF-binding sites or that the DNA-binding activity of ARF proteins might be promiscuous. In order to differentiate these possibilities, we tested binding of OsARF2 to the promoters of other sucrose transporters, including OsSUT2, OsSUT3, OsSUT4, and OsSUT5. Sequence analysis revealed that there are two SuRE elements in the OsSUT2 promoter, one in the OsSUT3 promoter, one in the OsSUT4 promoter, and two in the OsSUT5 promoter (SI Appendix, Fig. S9A). Yeast one-hybrid assay showed that OsARF2 could also bind to the promoters of OsSUT4 and OsSUT5, but not to others (SI Appendix, Fig. S9A). Consistent with these results, qRT-PCR experiments showed that the expression of OsSUT4 and OsSUT5 was also down-regulated in the Osarf2 mutant (SI Appendix, Fig. S9B). These observations suggest that in addition to OsSUT1, OsARF2 might also regulate the expression of OsSUT4 and OsSUT5 through direct binding to the SuRE elements in their promoters.
To further delineate the nucleotides in the SuRE motifs for specific binding of OsARF2, we conducted mutagenesis of each of the five nucleotides in the SuRE1 motif of the OsSUT1 promoter. Results showed that mutations of the first G, fourth T, or fifth C in the SuRE1 motif dramatically impaired the binding activity of OsARF2, suggesting that these three nucleotides are essential for specific binding of OsARF2 (SI Appendix, Fig.  S8 D and E). Interestingly, we also identified 8 SuRE and 10 AuxRE elements in the OsARF2 promoter (Fig. 4A), raising the possibility that OsARF2 might also regulate its own expression. EMSA experiments showed that OsARF2 cannot bind to the AuxRE elements but directly binds to the SuRE2 and SuRE7 elements in its own promoter, and the binding activity was gradually reduced by increased concentrations of unlabeled probes (SI Appendix, Fig. S10 A and B). Furthermore, luciferase reporter gene assay showed that OsARF2 could positively   , Fig. S10C). These results showed that OsARF2 activates the expression of OsSUT1 and itself by direct binding to the SuRE motifs in their promoters.
OsARF18 encodes a putative transcriptional repressor (37). A previous study reported that expression of OsARF2 was decreased in OsARF18 overexpression lines and that the OsARF18 overexpressors were defective in reproductive development, including abnormal flower and seed development (43). qRT-PCR analysis showed that OsARF18 expression was significantly up-regulated in the dao mutants, OsYUC4 overexpressors, and WT panicles treated with exogenous auxin (SI Appendix, Fig.  S11 A-C), suggesting that the expression of OsARF18 was induced by high levels of auxin. qRT-PCR and the histochemical staining of OsARF18::GUS (β-glucuronidase) reporter lines showed that the expression of OsARF18 was mainly expressed in stems, young anthers, and fertilized ovaries (SI Appendix, Fig.  S11 D-I). To further test the role of OsARF18 in regulating reproductive development, we generated OsARF18 overexpression lines in the WT background. We did not observe significant differences between the WT and OsARF18-OE plants during the vegetative growth period, except that the OsARF18-OE plants were shorter and had fewer tillers (SI Appendix, Fig. S12A). At the heading stage, the OsARF18-OE plants developed closed glumes and indehiscent anthers and produced parthenocarpic seeds (SI Appendix, Fig. S12 B-H).
To explore the relationship between OsARF18 with OsARF2 and OsSUT1 in regulating rice reproductive development, we compared sugar contents in the source and sink tissues of WT, OsARF18-OE, and Osarf2 mutant plants. Compared to WT, the sucrose levels in the sink tissues of OsARF18-OE and Osarf2 were all significantly decreased; however, higher levels of sucrose accumulated in the source tissues in the OsARF18-OE and Osarf2 mutant plants (SI Appendix, Fig. S13 A-D). Moreover, qRT-PCR analysis showed that expression of both OsARF2 and OsSUT1 was significantly decreased in the OsARF18-OE and Osarf2 mutant plants compared to that in the WT control (SI Appendix, Figs. S12H and S13 E and F). These observations suggest that OsARF18 might regulate rice reproductive development by repressing expression of OsARF2 and OsSUT1.
To test whether OsARF18 directly regulates OsARF2 expression, we performed transient expression experiments of a series of LUC reporter genes driven by different promoter fragments of OsARF2 (ProARF2-1 to ProARF2-4) in rice protoplasts. The LUC activity of ProARF2-1::LUC was significantly decreased, while the other three reporters (ProARF2-2::LUC, ProARF2-3::LUC, and ProARF2-4::LUC) had no significant change when coexpressed with OsARF18 (Fig. 4B). In addition, EMSA showed that OsARF18 directly bound to the EA1 fragment containing the AuxRE1 and AuxRE2 elements and this binding was gradually diminished by increased concentrations of unlabeled competitive probes (Fig. 4C). In addition, mutations of the AuxRE2 element abolished the binding activity of OsARF18, while mutations of AuxRE1 had no effect (Fig. 4C), suggesting that AuxRE2 is responsible for specific OsARF18 binding. Consistent with this, ChIP-PCR assay validated that OsARF18 could directly bind to various fragments of the OsARF2 promoter containing the AuxRE elements but not to fragments lacking the AuxRE motifs (Fig. 4D). Taken together, these results indicate that OsARF18 directly binds to the AuxRE2 element of the OsARF2 promoter and represses its expression.
Similarly, EMSA showed that OsARF18 could directly bind to the ES2 fragment containing the SuRE2 element in the OsSUT1 promoter and that the binding activity was gradually reduced by increased concentrations of unlabeled probes (Fig. 4E). Mutation or deletion of the SuRE2 element abolished the binding activity of OsARF18 (Fig. 4F), indicating that the binding was specific. As OsARF2 could also bind to the SuRE2 element to up-regulate the expression of OsSUT1 (Fig. 3M), we performed transient expression experiments in Nicotiana benthamiana to test the combined effect of OsARF2 and OsARF18 on the expression of OsSUT1. Coexpression of individual OsARF18 or OsARF2 with the P-S3::LUC reporter containing the SuRE2 element significantly decreased or increased the LUC activity of the reporter gene, respectively. When both OsARF18 and OsARF2 were coexpressed with the P-S3::LUC reporter, LUC activity was also significantly decreased (Fig. 4G). Taken together, these results suggest that OsARF18 could compete with OsARF2 for binding to the SuRE2 element in the OsSUT1 promoter, thus repressing OsSUT1 expression.
Finally, to substantiate the notion that reduced OsSUT1 expression contributes to reproductive defects in the dao and Osarf2 mutants, we overexpressed OsSUT1 in the dao and Osarf2 mutant backgrounds, respectively. Compared with these mutants, the spikelet opening and seed setting rates of the overexpressed transgenic plants were both significantly increased (Fig. 5 A and B, and SI Appendix, Fig. S14). Together, these results suggest that auxin plays an important role in regulating carbohydrate partitioning during rice reproductive development by regulating OsSUT1 expression.

Discussion
Understanding the molecular mechanisms regulating sourcesink carbohydrate partitioning is a fundamental question in plant biology, which also bears important implications for increasing crop yields (3)(4)(5). Although it has long been recognized that long-range transport of auxin and sugar from the source to sink tissues is important for growth and developmental processes, such as lateral root development, hypocotyl elongation, and shoot branching, little is known about the interplay between auxin and sugar (5,34,44). In this study, we uncovered an auxin-signaling cascade regulating carbohydrate partitioning during rice reproductive development. Several lines of evidence were collected. First, we showed that the dao mutant, OsYUC4 overexpressors, or WT plants treated with exogenous auxin are all defective in sucrose transport from the source tissues to the sink tissues, resulting in overaccumulation of sucrose in the source tissues (flag leaves) and stems and concomitant reduction of sucrose content in the sink tissues (lodicules, anthers, and ovaries) (Fig. 1 H-L). Second, we showed that the expression levels of OsSUT1 are significantly reduced in the dao mutant, OsYUC4 overexpressors, and WT plants treated with exogenous auxin (Fig. 2 D-F). Third, we showed that knockout mutants of Ossut1, Osarf2, and OsARF18-OE plants display a dao-like mutant phenotype and that they are defective in sucrose transport between the source and sink tissues during rice reproductive development (Fig. 2 G and H and SI Appendix,  Fig. S6 G-L and S12). Fourth, we showed that OsARF2 directly activates the expression of OsSUT1 and itself via binding to distinct SuRE elements in the OsSUT1 and its own promoter (Fig. 3 I-M and SI Appendix, Fig. S10). Fifth, we showed that OsARF18 directly represses the expression of OsARF2 and OsSUT1 by binding a highly similar AuxRE2 element (TGTC) and SuRE2 element (GTCTC) in their promoters, respectively  t test). (E) EMSA shows that OsARF18 binds to the ES2 fragment of the OsSUT1 promoter. (F) EMSA shows that OsARF18 binds to the WT ES2 fragments but not the fragments with mutated SuRE2. The SuRE2 elements are labeled in red. (G) Effects of OsARF18 and OsARF2 coexpression on luciferase expression driven by OsSUT1 promoter fragments P-S3 containing the SuRE2 element. Values are means ± SD (n = 3 biological replicates). Significance analysis was conducted with Student's t tests (**P < 0.01, ***P < 0.001, based on Student's t test).
( Fig. 4 and SI Appendix, Figs. S12H and S13E). Sixth, we showed that overexpression of OsSUT1 in the dao and Osarf2 mutant backgrounds significantly improves spikelet opening and seed-setting defects of these mutants (Fig. 5 A and B). Based on these findings, we propose an OsARF18-OsARF2-OsSUT1mediated auxin signaling cascade regulating carbohydrate partitioning between the source and sink tissues in rice, which is essential for proper development of rice reproductive organs, including lodicules, anthers, and seeds (Fig. 5C). Thus, our findings represent a major step forward in elucidating the regulatory relationship between auxin and sugar in regulating plant growth and development.
It should be noted that previous studies have reported that flowers of the Arabidopsis arf6 arf8 double-null mutants also fail to open due to lack of rapid swelling and enlargement of the petals at anthesis. These mutants also have undehisced anthers that do not release pollens (45). These phenotypes are remarkably similar to the closed glumes and undehisced anthers phenotype exhibited in the dao mutant, OsYUC4 overexpressors, and Osarf2 mutants. It has been proposed that the lodicules in grass plants are morphologically equivalent to petals in dicot plants and play an important role in flower opening (28)(29)(30). Thus, we speculate that a similar auxin signaling cascade may operate in both dicot and monocot plants to regulate flower opening and anther dehiscence.
Our results may also have important implications for breeding crops with increased yields. Several recent studies have demonstrated that manipulation of source activity or sink strength could significantly enhance crop productivity (4,21). For example, overexpression of the transcription factor Higher Yield Rice (HYR) associated with photosynthetic carbon metabolism promotes efficiencies in CO 2 assimilation and photosynthesis in leaves and improves grain yield in rice (46). Engineering of synthetic glycolate metabolism pathways or new chloroplastic photorespiratory bypass (GOC bypass: Glycolate oxidase, Oxalate oxidase and Catalase) increases photosynthetic efficiency and yield potentials in the field (47,48). However, the engineered GOC plants still suffer a low seed-setting rate. Thus, a major bottleneck to realizing grain yield increase is successful partitioning of photosynthates from the source tissues into the sink tissues. Our findings that auxin regulates carbohydrate partitioning during rice reproductive development may offer a viable approach to overcome this problem by better coordinating source and sink activities in crops.

Methods
Plant Materials and Growth Conditions. The background of various transgenic plants (OsYUC4-OE and OsARF18-OE) and CRISPR-Cas9 knockout plants (Osarf2-1, Osarf2-2, and Ossut1) used in this study is Oryza sativa subsp. japonica cv. Zhonghua 11 (ZH11). The background of dao mutant is O. sativa subsp. japonica cv. Nipponbare. All plants were grown in the Tu Qiao Experiment Station of Nanjing Agricultural University, Nanjing, China. All materials were planted at a spacing of 16.5 cm × 16.5 cm. A wide-row spacing of 23.5 cm was set between the plots.
In Vitro Pollen Germination. Pollen grains of WT and mutants were germinated on 1% agar medium containing 15% (wt/vol) sucrose and 20 mg/L K 2 B 4 O 7 . Spikelets were sampled just after anthesis and gently shaken above germination medium on a slide to collect pollen grains. The slide was placed in an incubator for 30 min at 28°C. Then, the number of germinated pollen grains was counted under a microscope.
GC-MS Soluble Sugar Assays. Soluble sugar levels were analyzed essentially as previously described (17,49). Materials including fresh anthers, spikelets, lemma/paleas, flag leaves, or stems were immediately frozen in liquid nitrogen and ground into powder. Six hundred microliters of methanol was added to the powder and vortexed briefly, then 50 μL ribitol (0.2 ng/mL, as an internal standard; Sigma-Aldrich) was added. The extracts were centrifuged at 10,000×g for 2 min, and the supernatant was immediately transferred to a new centrifuge tube and dried under vacuum for sugar assays. Acetylated derivatives of the sugars were made by mixing dry weight powder of plant sample in Me 2 SO with 150 μL acetic anhydride and 30 μL 1-methylimidazole and stirring for 10 min in glass tubes. Six hundred microliters of double-distilled H 2 O was added to the tubes to remove the excess acetic anhydride. Then, 100 μL CH 2 Cl 2 was added to extract the acetylated derivatives. The tubes were centrifuged for 1 min to partition the organic phases. Finally, 1 μL of the lower methylene chloride layer was injected for GC-MS analysis.
A GC-MS system (Trace-GC Ultra) connected to a Trace-DSQ mass selective detector with electrospray ionization set at full scan and selected ion monitoring (1.0 mass unit) (Thermo-Finnigan) was used to analyze the acetyl-derivatized sugars. The column used was a DB-17MS fused silica capillary column (30 m × 0.25 mm, 0.25-μm film thickness; Agilent). The GC oven temperature was programmed as follows: initiation at 100°C, then gradually ramping to 190°C (12°C/min) and holding for 6 min, subsequent ramping to 250°C (30°C/min) and holding for 6 min, and finally ramping to 280°C (40°C/min) and holding for 10 min. The injector and detector port temperatures were set at 250°C and 260°C, respectively. Helium was used as the carrier gas at a constant flow rate of 1.0 mL/min. An Xcalibur 2.0 workstation was used for data acquisition and quantification. 5, 6-CFDA Feeding Experiments. To examine sucrose uptake in the panicles of WT and mutant plants, feeding experiments using 5, 6-CFDA, a low-molecularweight dye, were carried out. The WT and mutant panicles were placed in an aqueous solution containing 5, 6-CFDA (100 μg/μL). After 30 min in the dark, the WT and mutant panicles were examined under fluorescence microscopy.
Radiolabeling Experiment. Sucrose radiolabeling was performed essentially as described (17,50). The WT and mutant flag leaves were placed in water containing 1 μCi [ 14 C] sucrose (21.8 GBq mmol-1 in ethanol:water [9:1, vol/vol]; MP Biochemicals). The stems were kept in this solution for 24 h at room temperature, using the stem submerged in water without 14 C-labeled sucrose as the control. Then, the anthers, ovaries, lemmas, paleas, flag leaves, and stems were collected for analysis. These materials were cut into pieces and incubated in 3 mL scintillation fluid. Radioactivity (counts per minute) was measured by liquid scintillation counting using a Beckman LS650. At least three biological replicates were measured.
RNA Isolation and qRT-PCR Analysis. Total RNAs were extracted from lodicules, anthers, ovaries, and flag leaves using TRIzol reagent (Invitrogen). Complementary DNAs were generated by reverse-transcription according to the manufacturer's instructions (Invitrogen). qRT-PCR was carried out using genespecific primers and SYBR Premix EX-Taq on a real-time PCR 7500 system (Applied Biosystems). Data were collected using the ABI Prism 7500 sequence detection system following the manufacturer's instructions. The rice Ubiquitin gene (LOC_Os03g13170) was used as the control. At least three biological replicates and three technical repeats were conducted. All primers for qRT-PCR can be found in SI Appendix, Table S1.
RNA-seq Analysis. Before flowering, total RNAs were extracted from the lodicules of WT and dao mutant using TRIzol and were purified using QIAGEN RNeasy Mini kits. RNA-seq libraries were prepared from WT and dao mutant with three replicates. The libraries were sequenced using Illumina HiSeq X-Ten platform (Biomarker Biotechnology Co.), and generated reads were cleaned and then mapped to the reference genome (Nipponbare). After normalization, gene expression levels were estimated by fragments per kilobase of transcript per million fragments mapped (FPKM). The multiple-testing P value <0.01 and fold change >2 were used to determine whether the gene was significantly differentially expressed or not.
Complementation of the Osarf2-1 Mutant. For functional complementation, an ∼6.8-kb fragment of genomic DNA containing the OsARF2 promoter region and the entire coding region was subcloned into the binary vector pCAM-BIA1305 carrying a hygromycin resistance marker to generate the p1305-OsARF2 construct. We induced homogenous Osarf2-1 calli, which were then transformed by cocultivation with Agrobacterium tumefaciens strain EHA105 carrying the p1305-OsARF2 plasmid and the control plasmid p1305.
OsARF18::GUS Reporter Gene Construct and Histochemical Staining. A 2.5-kb genomic fragment containing the promoter upstream of the ATG start codon of OsARF18 was amplified by PCR (primer sequences listed in SI Appendix, Table S1) using Nipponbare genomic DNA as the template and was cloned into the binary vector pCAMBIA1305 to drive expression of the GUS reporter gene. Transgenic plants were generated as described above. Images were captured using Leica Application Suite 3.3 and merged and enhanced using Photoshop CS (Adobe).
Dual-Luciferase Reporter Assay. To construct the effector plasmid, the fulllength coding sequences of OsARF18 and OsARF2 were inserted into the vector pCAMBIA1305.1-green fluorescent protein (GFP). For the reporter construct, a 1,800-bp upstream fragment of OsARF2 and a 920-bp upstream fragment of OsSUT1 were inserted into pGreenII0800-LUC to drive the firefly luciferase (LUC) gene to generate the ARF2 pro:LUC and SUT1pro:LUC constructs, respectively. These vectors were individually transformed into the A. tumefaciens strain EHA105. For transient coexpression, the effector and reporter constructs were coinfiltrated into the leaves of 4-week-old N. benthamiana plants for 2 d. Firefly LUC and REN activities were surveyed with a Dual-Luciferase Reporter Assay Kit (Promega), and the LUC activity, normalized to REN activity, was determined. The Renilla luciferase (REN) gene driven by 35S promoter was used as a normalization control. Transient transactivation with the reporter and the empty vector pCAMBIA1305.1-GFP was used as a negative control. Three biological replicates were performed for each assay.
Purification of Recombinant Protein and EMSAs. EMSA was performed using a light-shift chemiluminescent EMSA Kit (Pierce, 20148) following the manufacturer's protocol. The N-termini of OsARF18 and OsARF2 recombinant protein were fused in frame with glutathione S-transferase (GST) and expressed in the BL21 Escherichia coli strain. The fused proteins were induced with 0.5 mM isopropyl B-D-1-thiogalactopyranoside and incubated at 16°C for 28 h. The recombinant proteins were purified by GST-agarose affinity. The fragments of the OsARF2 and OsSUT1 promoter were amplified by PCR and labeled with and without 5 0 -biotin, respectively. The unlabeled fragments were used as competitors, and CKS-GST protein was used as a negative control.
ChIP Assay. The ChIP-qPCR assay was performed as described previously (51). Transgenic lines overexpressing OsARF18-Flag and OsARF2-Flag were used in this assay. ChIP was performed on flag leaves at the heading stage. The harvested samples were ground in liquid nitrogen. The powder was resuspended in buffer (50 mM N-2-hydroxyethylpiperazine-N 0 -2-ethanesulfonic acid [Hepes], pH 7.5, 150 mM KCl, 5 mM MgCl 2 , and 0.1% Triton X-100) supplemented with 1% formaldehyde and incubated for 10 min at 4°C. The powder was then incubated with 0.15 M glycine for 5 min at 4°C to quench the formaldehyde. The chromatin was subsequently isolated and sonicated to produce DNA fragments of around 300 bp. A flag tag-specific monoclonal antibody (Sigma-Aldrich, F1804; 1:300 dilution) was used for ChIP analysis. WT plants treated in the same way served as controls. The ChIP DNA products were analyzed by PCR using primers that were synthesized to amplify about 150-bp DNA fragments in the promoter region of OsARF2 or OsSUT1. All PCR experiments were performed under the following conditions: 95°C for 5 min, 40 cycles of 95°C for 10 s, 60°C for 25 s, and 72°C for 20 s. All primer sequences used for ChIP assays are listed in SI Appendix , Table S1. The experiment was repeated three times.
Data, Materials, and Software Availability. All study data are included in the article and/or supporting information.